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•IliQSis entitled s AlJIiAPlIONIC OSCILLA.COP.d IIT QUillTim ImJCHAHICS , 
oulmitted by SUBlDT' P, BiAlinAG-AP to the department of Physics, 
Indian Institute of xeclinology Kanpur, 
in partial fulfilment of the requirement of the Ph.h. degree. 


In this work the following eigenvalue problems have been 
inve s t i ga t e d j 

(i) tlie auartic erharmonic oscilla/tor (Hamiltonia.n ll=p -i-n h-Xx'', 


P = -1- 


d 


dx’ 


X >0) and the associated problem of fchc pure 


quartic oscillatoi- (H = p'^ + Xr'' , X > 0) 


2 , „2 , 


(ii) the general anharmonic oscillators (li = p t x + ^x 

^ >0, p = 3,4, . . . ) 5 

2 2 ^ 

(iii) the double minimuiQ oscillator (H p" •- x -i- Xx"’, X> 0) 


Eigenva.lues and eigenfunctions of these systems are 
obtained in all regimes of the quentum number n and the 
anharmonicity constant X , 'Ohe computed eigenfunctions are 
tlion used to obtain the transition momonts. the eigenvalues 
reported in this week eu-e accurate to 13 significant figures 
and the tro,nLition moments to 12 figax''OG. 

the eigenv.iluGS of the anharmonic oscillators 
(H = p^ + x^ -1- Xx^^'’') fall into two distinct classes (llioe 
Gt al. 1976) . In the low n, low X rogime the eigenvalues 
differ slightly from the harmonic oscillator levels whereas 
in the high n, hig^h X regime they differ slightly from the 



2 2 5 ' 

pure aniLarmonic oscillabor (H = oigenTaluoSo 

Betiroeii these t'^ro regimes lies the ‘h- .-undaiy layer' in which 
the eigenvalues ate neit]ier 'near harmonic' norc near 'pure 


anlia-rmonic ' . 


'h 0 existence of different 


re a ime s imp 1 i e 


different oscillation properties of the corresponding eigen- 
functions. Vje assert that t iis fa.ct must ho explicitly 
included in solving the eigenvalue prohlcra. '.Clio method 
applied in this work (Banorjeo .1976) involves the use of an 
appropriately (a.ccording to regime) scaled bat, is for the 
expansion of each eigenfunction. flie appropriately scaled 
basis simulates the different oscilD.ation properties of the 
eiA.onfunctions in different regimes and ma,kcs possible a 
uniform treatment of the problem in all regimes. 


f'ho Chapter I of the thesis is a review of various 
earlier methods used to solve the anharmonic oscillator 
eigenvalue px-ohlem. i'hesc methods arc suitable only in a 
pa,rticular regime of (n,X) a,nd do not gfve oigenvalacs to the 
same d,ccu37acy when oxtended to other regimes. '.Che construc- 
tion of an appropriatcl 3 ’' ocaJ-ed basis o,nd the metiiod used 
in fills thesis 'for the computation of the eig^envalueo is 
described in Chapter II. fhe actual computation of the 
eigenvalues is reduced to the determination of the roots of 
a transcendental equation in the energy. 'fhis is done 
numerically. Accurate eigenvalues and eigenfunctions of 
the quartic anharmonic and the pure quartic oscillator arc 



X 


thou o'bta.inod for vr.rioiis volucs of (iijX) covcrin/; a3.1 dif- 
foront regimes . In Cliaptor III >ro shoi: tliat the cigGiiTalue 
problem of the general anliarnonic or-icillet ors may bo 3o3.vod 
for cuiy ;i in cj.l rogimos of (n^X) iisinf_ cbe 3a,mo method. 

In Chapter IV accurate eigenvalues and eigenfunctions of 
the donble iiiinimum oscillator arc ca-lculatod. A IJEB expres- 
sion for the splitting botvrocn the lot/or eigenval\ies 5 bunched, 
in pairs, is obtained and the VJKB values a.re compa/od with 
the corresponding eiccurate values. In Chaptoi- Y, the tr'ansi- 
tion moments which arc the matrix olomorits of x^ (k = integer) 
between tjic anharmonic oscillator c;igens tartes are calculertcd 
usin£=' the computed eigenvalues :.nd eigenfunctions. Iho 
transition moments for any particvln.r transition satisfy an 
exciut linear recurrence relation (Banc;r3cc 1977) from which 
the higher moments for that transition may bo obtained 
recursively, without integration. 


Banorjeo K 1976 letters in Katli. Phys. 1 323 
Bane.rjee K 1977 Phys, Letters A63 223 

ITioo P I, IlacMillon D and IIontro.ll L ¥ 1976 J. Math, Phys. 

17 1320 



CHA.PTBR I 


iira?RODUco?iON 


Tlie study of anharmonic oscillators is a suliect of 
very genera3. interest. Tlie g^uartic anharmonic oscil3.ator 
desci’ibec. by the IIami3.tonian H = p + x +Xs;''^X>0 is the 
anharmonicity constant) is the simpD-est examjp^^e of nonlinear 
oscillation in c.lassica], mechanics (Krylov and Bo£;,oliubov 
3-945) . The invesxiga-tion of the same prob3-em. in quantum 
mechanics is usefu3. as a mode3. problem in mo3.ecTi3.ar physics 
and fie.ld theor 3 r In mo3.eculax vibrations, the potentiaD. 
functions are quite often of the mized harmonic-quart ic type. 
Indeed, in some cases of interest, dsie to accidentaJ. cancel- 
3.ation of the qua.dratic terms, the potentio..! functions become 
near3.y pui'-e quartic (see Cha.u a.nd Stelmeua .1963, Reid 1970). 

The quartic anharmonic oscillator is of particular interest 
in field theory because it is a. modeH. of quantum fie.ld 

olieory in one-dimensional space-time. The investigations from 
the fie-3.d theoxy point of viei; are referj-^ed to in the work of 
Bender and ¥u (3.976). A genera3. account on the anharm.onic 
OGCil-lator prob.lem may be found in Hioe and Montroll (1975) 
and Hioe, MacMillen and Montroll (3.976). 
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'Xliere have heen persistent attempts for finding the 
eigenva3,ues of the quartic anharmonic oscillator since the 
beginning of qn.antuin meclianics. In t}\c earliest attempts 
the formal- Rayleigh Schr‘c3ding<er pe.rbiirbation method ras 
applied to solve this eigenvaJme problem in which the Xx^ 
term was considered as the pertr.rbation on the harmonic 

p p 

oscil3.ator Hamiltonian Hq = p + Bender a,nd \m (1969) 

calcrhatcd nearly 70 expansion coefficients of the pert-orba- 
tion series in poT.roi's of X ioT the ground stare energy 
and observed an unexpected ra.pid increase in bheii magni- 
tudes. A. detailed analysis then revea,led tha.t the perturba- 
tion series for the eigenvaJ.ues in powers of X ±s not 
convergent for any positive value of X , no matter how 
small. Bender and I'iu (1969) studied the analytic properties 
of the eigenva-lues Ej^(X) (n being the quantum number) in 
the complex X plane and showed that there is a third order 
branch point at X = 0. It xras furtlier observed that X = 0 
is not an isolated singularity but is a limit point of an 
infinite number of branch points of E^(X), Bimon (.1970) 
proved the above properties of Bj^( X) bj^^ more rigorous and 
technical^.y sound arguments. In the language of perturbatioi 
theory the nonanalyticity of E^CX) implies that the perturba- 
tion series of B^(X) in powers of X is non- convergent. 

The following scaling argument due to iymanzik is 
si/nificant in this connection. On performing the scaling 
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• -1 

c r an s f o rma t i o n p-^a p on the Ilarailt onian H(ksA) 

2 2 

= p -( hx +Xx we obtain 

n(k,X) = a^x) . (l.l) 

a. 

Since t]ic a.bovo rjcalin^, bransf ormation io impleiuontahle by 

tinitary transf onmo/tion .tlio tiro aides of cnn, (l,.l) hove 

idcnticaJ. eij;],envoJ.UGS , 'Tlio relation ( 1 . 1 ) tl’us en:iti_e 3 

the the ei^.envalue prohlen of the qnaa'tic sniharnonic oscil- 

la!: Ox may be corn plot 0 I 37 doGcx“ibed in tome of blto rcducod 

2 2/ G 

llaiiiltonian H(lsX) = p + x -i- Xx h Sebtinp a'' = 1/X. one 
obtains 

syio) = . (i-2a) 

Hence , 

E^(l,l) 'v X-''/'-^ EjO,l) arj x-^co . _ (1.2b) 

3/3 

Since E^(0^3.) is indeponclont of X^, E^^(l 5 X) X"^ ^ for 

lar^o A . 

the difficulties in perturba cion approach may also be 
seen if one considers the equation ILfp = in momentum 

repre sento.tion (Hioe and Ilontroll .1975)° 

j 4 j 2 p 

( + p^) = E . (1.3) , 

dp ' dp 

The perturbation parameter X appears hero as the coeffi- 
cient of the hiyhest derivative . It is well hnovn that 
the expa.noion of the solnitions of such dif ferentia.l equsnbion 
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in power series of the small parameter is non-convergent 
(^anSyke 1964) . It is theref oi-e not so siirprising that the 
Rayleigh Dchr'6ding’er perturbation expansion of the anhar- 
monic oscillator eigenvalues in powers of X fai3.s to con- 
verge, Ihe perturbation treatment of this problem in clas- 
sical mecha,nics is instructive in this connection. It 
consists of obtaining approximate solutiors with toe help 
of the expansion of the displacement in power series of 
the anharmonicity constant. It leads to solutions that • 
contain sec^ilar terms like t^ sin at, t^ cos at in which 
the time 't' appears outside the sine and cosine symbols. 
Ihe sectilar terms in a finite term expansion introduce non- 
periodic solutions and cause the calcuJ.ated displacements 
to become arbitrarily large at 3.arge t. Further, the total 
energy becomes a function of time, violating the energy 
conservation principle (Bogoliubov and Iiitropolsky 1961) . 
fhe difficulties arising in the quantum mechanical perturba' 
tion expansion may be viewed in this context. 

A detailed analysis of the perturbation series of 
E^(X) shows that it is asymptotic in nature (Simon 1970). 
Such series are usually summed unique3.y through various 
summability techniques such as Stioltjes Fade or Borel 
methods. Loeffel et al, (1969) have proved that the 
perturbation series s\xms under Fade approximation to the 
actua.l eigenvalue. !IIhe Fade approximation in genera]. 
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consists of replacing the power series by a sequence of 
rational functions f(M,N) of form of a polynomia3. of 
degree M divided by another po3,ynomia.l of degree h. Simon 
(1970) ca3,culated the ground state energ37- by this approxima- 
tion for various va.lues of X . His resu3.ts sho'^v'- that the 
Pad4 approximants converge quickly for 0 < X <1 but for 
X >3. the rate of convergence is not very good. Graff i 
et a3.o (1970) described how improved va3.ues of E^(x) oan be 
obtained by using Pade approximants to the Horel summability 
method. Heid (1967) showed tha.t the perturbation series can 
also be summed bAr converting it into eouivalent continued 
fraction and obtained the gro\md state as we3J. as a few 
excited state eigenva3.ues lor various values of X. I'he 
agreement of the eigenvalues obtained from these various 
summabi3-ity methods with the corresponding accura,te eigen- 
values is found to be poor un3-ess n and X are suff tcj.entl5" 
smal3.. The scope of these methods in therefore 3.imifced to 
sma3.3. n, small x values only. 

Various variatiofial and numerical methods have been 
widelAJ' emp3.07ed by many authors either to the quartic 
enharmonic oscillator prob3.em or to the associated problem 
of the p-ure quartic oscillator (H = p + xx X > 0) , 
Ca3.culations are generally done in the harmonic o,scillator 
basis {x“ e~ ' } and the usual technique is t<fs truncate 
and diagonalize a large but finite matrix. Results of such 
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calcuH-ations are summarised in. iable (1.3.), Tliey differ 
from one s-notlier either in the method of diagona.liza-tion or 
in the size of the basis used. Some cii-bculations for upper 
and D-Ower bounds of first few eigenvalues with varying 
anharraonicity xrere done by Bazley and Pox (3.961) and .Reid 
(3-965). fhe procedure for 3.ower bounds used by Bazley 
and Pox was to construct intermediate llami-ltonians such 
that ... <H and to determine eic-cnva.3.ues of 

successive H . Reid used the method of L'dwdin (1965) for 
obtaining the lower bounds. The upper bounds were ca3.culatod 
by emp.loying the usua,l Ilayleigh-Ritz vai-iational approach, 

It may be noted that in these calculations the gap between 
the upper and the lower bounds increases rapid3.y on increasing 
n or X (e.g. in Reid's work the gap, which is 0(10“'' ^) for i 
B(-( X = 0.25), becomes 0(.10“^) for E^( X = 1.0)). I 


Biswas et al. (1975) used the 'Hi3..1 determinant' method | 
to obtain first eight eigenvalues of tlie quartic anharmonic | 
oscil.lator ior values of X in the rang'e 0< X^IOO. They 1 

2 i 

used an expansion in terms of the functions { x'^ e“^ I 

for the eigenfunctions and obtained eigenva.l.ues by finding 
the roots of a sequence of characteristic polynomials in E. 

The polynomials corresponded to various order truncations 
of the Hil3. determinant. Por higher eigGnva3.uos or for 
higher X the numerica.l errors in their work become too severe 
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iriae results of the variational and numerical calcula- 
tions show that these methods are suitable for the evaluation 
of a few lower eigenvalues and for moderate anharmonicities , 
The accuracy of the results is seriously affected on increas- 
ing n or X and is not necessarily improved by increasing 
the size of the basis. Further, the eigenfunct ions obtained 
from variational calculations are far .less accurate Than the 
corresponding eigenvalues. Evaluation of matrix elements of 
operators using such eigenfunctions is not expected to yield 
values of known or definite accuracy. 

% 

The ¥EIB method has been used for obtaining approximate ; 
eigenvalues for high n. In the ¥KB approxima.tion (Titchmarsh 

1961) , ■ I 

1 

E^(0,1) = C + C =22/^(1.376) .(1-4) | 

I' 

It is known from eqn, (1.2b) that the quartic anharmonic I 

— 1 / 3 i 

oscillator eigenvalues E^^(3.,x) = \ ' E^(0,l) in the large j 

X limit. Hence, for large n, large X - | 

Ej^(l,X) = C X^/^ (n + 1)"^/^ . (1.5) 

A more sophist j cated ¥KB expression for E^(l,x) is obtained 
by Hioe and Montroll (1975); see eqn. (l.6.c) later. 

The most comprehensive work on the anharmonic oscilla- 
tors is due to Hioe and Montroll (.1975) and Hioe et al. (.1976 
They distinguished two limiting regimes of values of n and x- 
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In one regime the energy eigenvalues differ slightly from 
the harmonic oscilla.tor levels; in the othex they differ 
slightly from the pure quartic oscillator eigenvalues o '.Chese 
regimes are called the 'near harmonic' regime and the*near 
quartic' regime, respectively. Between these two regimes 
lies the 'Boundary layer' in which the eigenvalues are 
neither 'near harmonic' nor 'near quartic', Ihe^^- developed 
fast converging algorithms for computing the eigenva-lues 
in small n regime hy writing bhe eigenvpJ.ue prohlem in 
Bargmann representation and nolving, fho associated difference 
equa.tion. The eigenvalues were thus computed to 8-9 signi- 
ficant figures for n = 0, 1, , , , 8 and for values of X in 
the range ,004^ ^ :S40000. They also constructed several 
simple formulae for with different ranges of validity 

which when combined give good approximations (about 8 signi- 
ficant figures) to ) except in the 'boundary layer'. 

The formulae with their ranges of va.lidity are as follows; 


(a) Bor the near harmonic regime 


E (X) 
n^ ' 


(2n+l) + ^X { l+2n(n+l) J 


2 ( n+1 ) ( n+^) ^ ( n+ 2 ) 



[4+3 ^(2n+3)3 




32r4+TxT2n+^ 


n(n-^)^(n-l) 

T4+3x^r2n^Tl 


32[4+3X (2n-3) ] ? ^ ^ ^ ^ 


( 1, 6a) 
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(■b) Por small n, lojrge X region 

1^(X) [s^ + + p^X + ...], (1.61) 

-Jbere tlie constar.ts s , a and B were determined bj fitting 

n n '^n 

eqn. (D..6b) to the nimerical .values of 3^(X) for ea-clx stace 
up 1 0 n=3.0 . 


(c) Por leo.-’ge n^ large A region 


E^(X) =X-‘^/5 


[ c { ( n+i) + 


4/3 


+ b X' 


2' 
4/3 


( n-:-A-) 


} + a ( 11 “!“^) 


1^2/3 , -2/3 


+ . . . ] 


(l.uc) 


The above oxproosion reou.lts from a detailed investigation 
of the ¥E:b approximation formula and consists of expounding 
bhe o3.1iptic integrals in the ¥KB formula in a series in j 

the constants c, a- and b arc identified from tliis expansion. i 

It may be mentioned hero that Mathe’^s a.nd .'Jswaran (1972) alsoj 

[, 

ol'tainod some approximate formulae throug.h a semic3assical treatmai 

i 

Rocenbly, Iiav.l-io .'rs o.nci C-ovinda.r-ajan (1977) used a 'residue 
squaring method ' foi’ che iterative d i agonal i sat ion of the 
ciuouctic anharmonic oscillacor hamilt onian in which the Xx^ I 

s 

term is ass-uriod to be a pertarbation on the rest. In this i 

I 

work the off diagonal part (of order x compaxed to the | 

2 ^ 8 ^' 

diagonal part) is successively reduced to oi*c;Qra X , X‘, X » I 

i 

r 

...... They obtained foti;r lovrest even parity eig-envaluos and | 

|! 

an approximate ana.lytic formula for E^(x ) , similar to the | 
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eq_uation (1.6a), x^hich. £ives i^ood approximations to the 
eigenvalues in the nc-ar hai’monic regime. 

It is seen irom the above discrission bhat the various 
methods, XThicli have boon applied to solve tho anharmonic 
oscillator eigenvalue problem, are suitable only in some 
particular regime of (n,X) and do not give eigenvalues to 
the same accuracy uhen extended to other regimes. Moreover, 
th 03 do not yield accurate eigenfunctions. In contrast, we 
use a method (Banor.jee 1976) to obtain the eigenvcluos of 
the quartic anharmonic oscillator x^hich applies x-j-ith xmiform 
and arbitrarily high accuracy for all values of n and X . 

I'liG method also yields eigenfxmctions of accuracy comparable 
with that of the eigenvalues which are used for the compuba“ 
tion of high accuracy matrix elements. In the next chapter 
the method is described and the eigenvalues and the oigen- 
fxnctions of the quartic anharmonic oscillator and the asso- 
ciated problem of the pure quartic oscid-lator arc obtained. 

The eigenvalues, accujrato to 15 significant figures, are 
presented for various values of {n,X ) covering all tho 
different regimes. ¥e shox-r in Chapter III that the method 
may be extended to solve general anharmonic oscillator 
(H = p^ + x^ + Xx^*^, p, = 3,4s ...) eigenvalue problem. 
Eigenvalues of the sextic ( p=3) and the octic (ii=4) anharmonic 
oscillators are thus obtained accurate to 15 significant 
figures in all regimes of (n,X). In Chapter IT the eigenvalue 
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2 2 4 

spectrum of bhe double minimum oscillator (ll = p “X +Xx , 

X > 0) is investigated. 'I'Lc lo-wci' eigenvalues of the double 
minimum osci3.1ator are c3-Ose3.y bunched in pairs loi* small- X. 
iiiose eigenva.lues have been eva.liiated accurately using the 
samo method. A ¥KIB expression for the splitting between 
the eigenvalues bunched in pa,irs is obtained and the ¥KB 
values are compered with the corresponding accurate vcJ.uos. 
In Chapter V the transition moments between the anhexmonic 
oscil3.ator energy eigenstates are obtr^ined from the computed 
eigenvalues and eigenfunctions. Further, the multipole 
transition moments are shoim to satisfjr an exact 3.j.noar 
recurrence relation which is va.lid for a.ny po.lynomial 
potential. 



CHA.PTEE II 


'IHB QUAPtllC AITIIAEJ-IOFIC OSCILLATOR 


In thic' chapter ve determine accurate eigenvalue c and 
eigenfunctiona of the quartic anharmonic oscillator 


(H = 4- ^ ^ Oj p 


dx 


) for various values of 


the quantum number n and the anharmonicity constant The 

existence of two distinct regimes of values of (n>X)^ namely 
the 'near harmonic' and the ‘near quartic', separated by a 
'boundary layer' implies different oscillation properties of 
the corresponding eigenfunctions. ¥e assert that this fact 
must be explicitly included in solving the eigenvalue problem, 
following Banerjee (1976), it is shown in the next section 
that the basis functions (used in the solution of the eigen- 
value problem) may be appropriately sca3.ed to simulate the 
oscillation properties of the eigenfunctions in all regimes. 


II . 1 Scaling and t he Appropriate Sc a ling Formul a 

The eigenfunctions '{'^.(xjX) may be expanded in the basis 
2 ^ 


/ m„-ax -i 
IX e / 


as; 


\) 


-ax 


m=0 


m 

a X ^ 
m ^ 


( 2 . 1 ) 
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where the scaling is introduced through the pars^meter 'a'. 

For an eifecti-ve expansion the scaling c is cb.osen such that 
a sufficient nutaloer of t]iG lo^'er laemhers of the basis func- 
tions (at least n for the n-th state) have their main contri- 
bution in the region of oscillation and outside i'c thev decay 
iiionotonically Since tho region of oscillation depends on 
n and X , the scaling a also depends on n and ^ . A simple 
criterion for the appropriate scaling is obtained as follows: 
fhe region of OGCil3.ation for the n-th eigenfunction of the 
guartic anharmonic osci.1 la.tor (f or sufficiently' lange n) : 

In the ¥1IB approximation ) ~ CX^''^^(n + ^ , 

therefore the region of oscillation The • 

exact n-th eigenfunction has n zeros in the region of occil- ! 
lation. Hence setting the region of oscillation for the I 

n-th eigenfunction equal to the width of the n-th basis ' 

___ " I 

function {'^'({n/a)) we obtain | 

i 

! 

a(n,X) 'V* (n + X^/^. | 

li'iis puts tho span of the first n ( or a number pi-oportional 

f 

to n) basis functions in the region of oscillation of the 
n-th eigenfunction f o r al l n and X , just a,s required for an 

i 

effective expansion^ In view of the ¥KB estimate and the | 

large X assumption implicit in tho derivation of the hhove i 

scaling formula, it is not expected to be good when n and 
or ^ is small. However, for X -»■ 0 or for small n the 
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scaling must approach the va.lue 1/2 apijropriate for the 
harmonic osci.llator. Hence the Dca.ling formula for all 
regimes of (n^X) is 

a(n,X) = I -1- (n + X^^^ . (2.2) 


The effect of using an appropriate 13?- scc.led basis is remarkable. 
It is now possible to compute the eigenvalues in any regime 
of (ujx) with arbitrarily high accuracy. 


It is to be noted that the different regimes of (n,X) 
ai’e distinguished according to the a.bove sca]ing formula. 


Thus 


1)1/5 g/3 1 


<< 


2 tlie near liarmonic regime^ 


(n + i) 


lvl/5 ,1/3 


>> 


is the pure quartic regime. 


(„ + 1)1/3 g/3 . 1 


2 is the boundary la3?'er between 
the above two regimes. 


The value of the combination (n + X is seen to determine 
the regime to which an eigenvalue be.longD, The importance of 
the above combination of n and X in determining the various 
regimes was also recognised by Ilioe et al, (1976) on essen- 
tially empirical grounds. It is in this work that the 
combination (n + ^) X is shoim to determine the characteristic 
scaling in a given regime of (n>X) through the relation (2,2).- 
This observation leads to the construction of a sca.le adapted 
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basis and makes possible a unifOa.’m treatment of the problem 
in all regimes of (n,X), 

II. 2 Method 

fhe method applied for obtaining the eigenvalues is 
described in this section. 1‘he dchrbdingor equation fox’ the 
quartic anharmoni-c oscillator is 
2 

[- + Xx^’] i|;(x*,X) = B(X ) (2.3) 

dx"^ 

uhei'e the eigenfunctions -»■ 0 as x -+• 4^ . flie expan- 

sion (2.1) on substitution into the above equation yields the 

Collowing 4“temn linear recurrence relation, connecting the 
•alternate expansion coefficients > 

(m+l)(m+2) a ^_^2 + (E“4o:m-2a) a ^^_2 ^ ^ni-4 

( 2 : 4 ) 

■flic ebove recurrence relation may bo rewritten in the follov— 
ing notations 3 

®'m+2 ^m,m®'m , m“2®'m“2 ' ^m,m“^^m-4 ~ 

uhero 

^ = (E - ^4..am-2q:l , ^ 

^m,m ua+T) (m+2) * m,m“2 X^+lTT®?^ ’ 

^m,m“4 " (mflRnH-^ 


(2.5) 
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Since the Hamiltonian for the system has even symmetry, the 
solutions of the Sclartidinger eq_uation are either even or 
odd functions of x. fhe even and odd parity solutions are 
obtained respectively by assigning the initial conditions 
(i) a|-^=l 5 a-|=0 (ii) 8.q=0, a^=3.. fhe recursion( 2 . 5) may he 
vieued as an infinite set of linear homogeneous equations 
in the unlcnoT/ms {a }. lor the self consistency the infinite 
determinant &(E) formed from the coefficients of ^ must 
vanish. It gives 


A(E) = 


o 

o 

1 

0 

» • 





^20 

^^22 

1 

0 

« o 




^40 

<^42 

"^44 

1 

0 

« • 



0 

^^62 

"64 

, "66 

1 

0 



0 0 

0 

"m,m“^ 

d „ 

ni,in-2 

d 

m,m 

1 

0 

« • 

0 0 

• « 

• « 

m 0 

« c 

m <• 

0 B 

• 0 


= 0 .( 2 . 6 ) 


ihe roots of t]ie above transcendental equation are the eigen- 
val-ues. This mode of writing the char8,cteris tic equation in 
the form of an infinite determinant is wel.l knoxm from the 
eigenvalue problem associated with Hi3J. 's equation ( l/hit taker 
and Hatson 1927) . Denoting the trmicated determinant formed 
by omittii^ all rows and columns beyond the element ^ 
as s noted, that is a polynomiaD. in E 
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of degree m/2 when, m is even, and (m~l)/2 when m is odd. The 
zeros of niimerically stable for 

large m pi'^ovidod the scaling a io appropriate. The deterini- 
nant m-><» defines A(E). Hence the problem of 

ob 'Gaining the eigenvalues is reduced to finding the zeros 


ox Aj^(E) for sufficiently large 


m. 


The zeros of Ajj^(E) are 


obtained numerico.lly by Hewton's method which requires an 
initial estimate of the required zero say well 

as the valu.es of 'A^(E) and A '(E) (= 9A /3E) at E = E. 

To evaluate note that the truncated determinants 

satisfy the following 4“term linear recurrence relation 

A^(E) + A^_2(E) - A^_^(E) =0, 

(2.7) 

obtained by expanding the determinant veJ-ue of 

the determinants A^(E) ma.y thus be computed successively upto 
any order in terms of A^(or A^.^) using (2.7). The recurrence 
relation (2.7) on differentiation wi ch respect to E yields 


m+2 


(E) - d' A (E) - d (E) A 1(B) + d .. a’ o(E) 
‘ m.m m^ ‘ nomi ‘ m^ ^ m,m“2 ^-2^ ’ 


d 


m,m“ 




0 . 


( 2 . 8 ) 


from which Hiay be computed upto any order recursively. 

The recursions (2.7) and (2.8) are numerically stable. The 
initial estimates for the eigenvalues required in the 
Hewton's’ method may be obtained for low n by evaluating a 
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sufficient!’/ large order determinant Aj,-(E) from the recur- 
Sion (2.7) £it various E points. Opposite signs of Aj„(B) 
tor t'lvo neighbouring B valties indico,te that an eigenvalue 
is crossed which provides suff icieivLly accurate estimate for 
the Bewton's met/iod. Bor high n (and' not too low X) the 
coiresponding approuimationo of the eigenvalues are good 
initial estimates. 'die procedure fo3r obtaining' initial 
estimates in the (high nj, low X) region of the 'boundary 
3.8.yer' is described in Appendix; A. 

The actual computation of the eigenvalues may now be 
performed in fche following manner. An initial estimate 
(B = ^j_ 23 ^j_-ti 8 l^ into the recursions (2.7) and (2.8) 

containing the appropria,te i^lue for a. The recursions are 
then continued on a computer until the corrections ^ given by 
the Newton's formula 

\ 

6E(m) = -{ A„(E)/a;(B)} _ , (2.9) 

^--^initial 

stabilize to a prescribed extent ( see Wilkinson 1965) . The 

corrected value fo?? B is then fed back in the second step 

as the initial value. This is continued till the required 

accuracy in the computed eigenvalues is reached. Due to 

the quadratic converg’ence of Newton's method it is possible 
* 

to refine a rather crude initial estimate for an eigenvalue 
(say, within a few percent) to a 15 figure accuracy in 4 or 
5 steps for all eigenvalues. Some typical examples to 
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elucidate this are given in Table (ll.l). In principle, 

Che method can be carried to a.n arbitrarjr high, accuracy. 

The accuracy of the computation is 3.imited only by the preci“ 
Sion of the arithmecic used (.16 significant figures in 
IBM 7044 ) . 

II . 3 E igenvalues 

Very accurate eigenva.lues of the auartic anharmonic 
osci.llator and the associated problem of the pure quartic 
oscillator (I! = p + Xx ‘ , X > 0) have been obtained using the 
method described in the previous section anci the sca.ling 
formula 

a(n,X ) = I + (1.2 to I. 4 ) (n + X^/^. (2.10) 

The constant xrithin the bracket (1.2 to .1.4) has been set j 

empirically by finding the values of the scaling a for which 
che computed eigenvalues stc!,bilize the earliest. The constant 
(.1.2 to 1 . 4 ) is found to work admirably in the entire range 
of (n^X). The ca.lculation of the eigenvalues ior the pure ! 

I 

quartic oscillator is similar to tha.t done for the quartic j 

I 

anharmonic oscillsfcor. The only difference is in the value | 

I 

of equations ( 2 . 5 ) to (2 . 8) , which in the pure quantic ^ 

' \ ^ 9 

oscillator case is equal to (4a /(nH-l) (m+2) ) . The computations 
were done on IM ’7044 computer using double precision arith- 
metic (16 digits mantissa). The eigenvalues were evaluated to 

^ Tables referred to in any chapter are given at the end of 
that chapter. 
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16 significant figures and then rounded off to 15 figures for 
the i'ableSo 

The eigenva3.ues of the quartic anliarmonic and the pure 
quartic oscillator for various values of (njX) covering all 
regimes are listed in 'Tables (II. 2 to II. 5). In Table (II. 2) 
ve compare om' results of E^(X =1) for ii=0j 10^ 100, 1000 and 
10000 with the corresponding resuJ.ts of various earlier 
calculations. In Table (II, 3) the first 50 eigenvalues of the 
quartic enharmonic and the pure quartic oscillators are given 
for X=l. The eigenvalues for any other value of X in the case 
of pure quartic oscillator can be obtained from the corres- 
ponding values for X =1 through the exact scalinjg relation E^( X) 

= X''-^E^(1). Eor the quartic enharmonic oscillator different 
eigenvalues for various values of X between .00001 and 40000 
iTere computed, and are presented in the Table (II. 4). It 
includes the eigenva.liles for the (high n, low X) region of 
the 'boundary l^ayer' . The computation of eigenvalues in this 
region of the boundary layer is found to be the most difficult 
in the earlier literature. Hone of the approximation formulae 
(1.6as,b,c) constructed by Hio^ et al. (1975) are adequate for 
this (high n, low X) region. In the present work the eigen- 
values in this region (e.g. eigenvalues corresponding to 
n=100, X slO”^ and n=1000, X=10*”^ in Table (II, 4)) are obtained 
by the same technique and with the same accuracy as any other. 
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liigh3-iglitiiig the scope of computations with an appropriately 
scaled "basis. In tafolo (II. 5) we focus on the regimes of 
extreme values of (n^X). Hioe et a3.. (1975s 1976) give 
different formulations of the eigenvalue proolem in each of 
these regimes wliich cannot be extended into other regimes 
because of the boundary layer in between. In contrasts, ve 

have obtained eigenvalues in each of those regimes by the 
same formu3.ation and with the same accuracy. 


II. 4 Stability of Zeros of Aj^(5) 


The recursive evaluation of the determinants A (e) and 

m 

the stabilit;?' of their zeros will noxr be considered in some 
detail. It may be noted that the recursion (2,7) is obtain- 
able from the recursion (2.5) by replacing a^^^ with A ^ and 
changing the sign of every alternate term. This px^escription 
is valid when the coefficient of the highest ore'er term in 
the { } recursion is set unity (by properly dividing, if 
ne ce ss ary ) , Then 


^m = (-1)“'^^ a^ , m = 0,2,4,..., 

= (-l)^'^"^^/2a^, ni=l,3,5,... . 


( 2 . 11 ) 


In the case when the coefficient of the highest order term 
is not unity, the (a^} recursion is 

■# 
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liigliD-igliting the scope of computations with an appropriately 
scaled basis. In ’.'able (II. 5) wo focus on the rogimes of 
extreme values of (n^X). Hioe et aJ.. (1975s 1976) give 
different formulations of the eigenvalue proolem in each of 
these regimes wliich cannot be extended into other regimes 
because of the boundary layer in between. In contrasty wo 

have obtained eigenvalues in each of those regimes by the 
sa,me formulation and with the same accuracy v 


II. 4 Stability of Zeros of 

The recursive evaluation of the determinants and 

the stability of their zeros will now be considered in some 
detail. It may be noted that the recursion (2.7) is obtain- 
able from the recursion (2.5) by replacing with A ^ and 
changing the sign of every alternate term. Thin prescription 
is valid when the coefficient of the highest orfer term in 
the { a^ } recursion is set unity (by properly dividing, if 
necessary) . Then 


^m “ (-1)^'^^ a^ , m — 0,2,4,. 

A = (~l)^“"^^/^a , m=l,3,5,... • 

■m ' ' m 


( 2 . 11 ) 


In the case when the coefficient of the highest order term 
is not unity, the recursion is 

# 
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HI y 

ni 4-2 

Ql J-H 

m +2 

+ d 

m. 

m^m '^m,m“ 2 ^m “2 

d .a . = 

m,m "4 21-4 

'file 

corresponding 

infinite determinant formed in ' 

io 



'^02 

0 




*^20 

^22 

^24 





^A 2 

^44 *^46 ^ 

# O 

A(S) 

=: 

0 

« « 

1 

%2 

« » 

^^66 *^68 

« • * 0 *0 

0 .. . . 

• 0 . « « « o 



• • 

0 

d . d o d ■ 

m ^-4 m , m -2 m , m 

^-m,m-i -2 ° '• 



• « 

* • 

« o • • • • 

* o # o «t c 


0 . (2.12) 


. (2.13) 


■Xl'ie trunco.ted determinants formed by omitting all 

rovs and coD.-amns be37'ond the element d^ in S(E) may be 

^ I 

expanded in terms of the tmnca'ced determinants of 3.ower ; 

orders. It provides the following 4-term recnrrence relation^' 

! 

I 

- <^m,m + V2,m'^m,m-2 I 


5 m--2 ^m-2 , m^m , ^ m~4 ^ ^ ^ ’ 


( 2 . 14)1 


The determinants {A^(B)} are related to {A^(E)} by 


\(E) - cljj^_2,m ‘^ra»4,m-2 *^02 ’ m-0,2,4,. 


= V2,m '^m-4,m-2 


( 2 . 15 ) 


« o o 


^ m—l , 3 y 5 7 * '» 
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The asymptotic behaviour of the solution of the 
Schrbdinger equation (2.3)'^ exp ( + lxj^/3). It requires 

%+2/'^m-4 X /m^ , m - - , ' (2.16) 

in the se3ries solution (2.1). It follous immediately from 
the relation (2.11) that for the doterminanto {A^ it is 
re quired that 

^m+2/\i-4 > m ~ . (2.17) 

IJe observe during the recursive evaluation of the value of 
^m+2^^^ that for sufficiently , large m the second and third 
terms in the recursion (2.7) become order of magnitudes ; 

smaller than the last term and the asymptotic dependence 
(2.17) is satisfied. The asymptotic relation (2.17) implies | 
a sequence of decreasing determinants beyond a sufficiently 
large m. This zeroing for large m must be isolated from the | 
determination of the eigenvalues which are the zeros of A(E) | 
for values of E. It is possible to achieve this quite I 

simply by redefining the determinants ^jj 3 ^(E) after multiplyin^j 

1 

them with a large number whenever, while applying Eewton's 

method, the recursively computed determinants become too 

small in magnitude. This renormalization amounts to starting 

the recursion with a higher value of the arbitrary constant 

Aq (or A^) . More generally, other recursively connected 

sequences of determinants {A„(E) J may be defined such that 

ni 
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the zeros of and are coirnion hnt may he 

f'iven any desired as 3 ?’mptotic behaviour for larye m. fliia 
is done hy multiplying the recursion (2,5) by a function of 
m oa}’- f(m). The corresponding infinite determinant a(B) 
and its various order truncations are related to the respec- 
tive quantities for a(E) by the relation 

IJE) = f(m-2) f(m-4) ... f (o or 1) (2.18) 


vhere tbo ihs contciins f(0) (or f(.l)) for che even (or odd) 
eigenvalues. Clearly can bo given any asymptotic 

behaviour for large m by properly choosing the fimction f(m). 
bince f(m) is independent of B b,y definition, the zeros of 
Aj^(E) and Ajj^(E) are common and they are equally suited for the 
computation of the eigenvalues, ‘the 'renormalization' 
suggested above is a special, case of this multiplication in 
vhich all rows are loft intact except one which is multi- 
p.lieci by a large number. ' 

To see the stability of the zeros of ^j^(E) as m-^ " 

''■TQ firsb consider the same problem for the harmonic oscil- j 

let or (n = p +x ) whieJi is exactly solLible. In the case of J 

! 

t} 3 e harmonic oscillator the zeros of the characteristic | 

i 

polyaomial A (E) are real and the ratio of successive i 

m I 

\ 

polynomials is ! 


V2(®) 


A (B) 

m'' 


(2m+l - E) 

feDWfy 
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■The zeros of coreiet of all zeros of \(E) plus a 

zero at E = 2nH-l. Thus the eigenvalues res'olting from the 
solution of A_^(E)=0 are reproduced ero.ctly 'by tlie solutions 

charsccteristic of an exactly 
solr.hle problem, for tne qu:%rtic anharmonic or the pure 
ouartic oscillator the zeros of successive order po3_ynomials 
are different. I'or a. given JJ and sufficiently large ni; 


however. 




^rji-4^ 3) 


m ' 


m 


mi 

i 


lib implies that tie largest order term in the ratio of 

independent of B for m-*°° . The successive 
polynomials as a function of E therefore differ by a multi- 
plicative constant ( depending upon m) for .large, m. Hence, 

the zeros of A^(E) stabilize for large m. 'phe above discussioni 
on stability folloi/s Eanerjee et al. (to be published). 

11.5 Checks for the Eigenvi.lues 


The confidence in the accuiacy of the computed eigen- 
va.lues is derived from the i allowing checks; 

(i) Computations were done with several different initial 
estimates for each eigenvalue. The intermediate numbers 
invo.lved in the computation are different for different 
initial estimates but the final results for 'che eigenvalues 
remain the same in all cases. 
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(ii) The eigenvalues were computed using severaJ. vaJ.ues of 
the scaling a in the range given bjr ecj_n. (2.10). Its effect 
is that tie stabilization of an inh-tic''! estimate for an 
etgenve.lue to 16 figures occurs Oit slight?-y different point 
in recursions (2.7) and (2,S). Each l'^(X) obtained in this 

0 01 * 1 : was thus checked and conlirmed by using 7 or 4- different 
values of k. 

(iii) Three separate computo.tions using an increasing!, a 

decreasing and a nearly flar sequence of determinants (see 
section II. 4) yie-lded the same valu.es for ) • 

(iv) bLifficiently large order determinant containing the 
.appropriate scaling is evaluated for two neighbouring values 
of E, Opposite sigiS of indicates that an eigenvalue is 
crossed. The computed eigenvalues are thus tested and in 
bhe process are upper and lower bounded in the last signi- 
ficant figure. 

II. 6 Eigenfunctions 

vJhen S is set equal to a computed eigenvalue in 
recursion (2.5) the resulting coefficients {%(S)} provide 
a very convenient representation for che corresponding 
eigenfunction through the expansion (2.1). The following 
properties of the solutions of the Schrddinger equation 
iKxjE) are important to note in this connection (Titchmarsh 


1961 ) ; 
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(i) Tbe eigenfunction associated x^itli the n-th eigenvalue 

has n jseros in the classicaD. region and 0 

as well as ) -^ 0 as x ■♦■-roo . 

11 iT 

(ii) lor E not exactly equeJ_ bo an eigenvalue,'!^ (x;E) can 
have at most one zero in blie nonclass ica.1 re:^,ior; on both 
sjdes of tl±e axj.s and as x -^ + <» it goes either to +«> or -«>. 

( iii) 'J.'iie blow up of \J>(x;S) starts in the nonclassical region 
and shifts to larger [xf as E approaches an eigenvalue. 

(iv) As S crosses an eigenvalue, ’!> (xjE) changes its sign 
of blow up in the nonclassical region and tends to infinity 
irith opposite sign. 

fhe expansion coefficients {a^^CE)} are evaluated 
recursively from (2.5) for various computed eigenvalues. 

It is observed that the sign of { a^} stabilizes after a 
certain sufficiently large index either to plus or minus 
sign w'lich implies that the compr.ted 4)(x,E) ->+oo or -<» as 
x-^ + oo. However, this large |x| heh.avio-ui' does not affect 
the computation of the accurate eigenfunctions significantly. 
Since the eigenvalues used are accurate to 15 sirpaif leant 
figures, the computed eigenfunctions reach extremely small 
values in the nonclassical region before the blow up starts. 
•Qho part of the computed eigenf inactions where the blow up 
occurs for large |x | may therefore be replaced by zero 
without losing much information. 
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'i'o test how well thu computed eigenfiirLctions satisfy 
bhe Schrddinger equation we compare the bwo sides of 

x) /i|^( x) = E at various points x„ For the first ten eigen” 
Tunctions which we have computed for X = i the test equality 
is satisiied to 19”1<- significant figures from x=0 to points 
■troll outside the classical region. For example, in the case 
of tie loth eigenfunction of the quartic anlaarmonic oscilla” 
tor with Xr.1^ the test equality q/ ''^ 0 “ ^10 satisfied 
to at least 19 significant figures in the entire classical 
region (x At a point x = 1.5 the test 

equality is still satisfied to 10 significant figures, where 
the value of the computed eigenfunction (x - 1.5 

ir u(10“^^) relative to i|'_^q(x= 0)=1. flie accuracy of the 
computed eigenfunctions are also checked by evaluating them 
"i.'or two neighbouring; values of E which upper and lower bound 
the eigenvalues in the 15th signi^'icanx figure. Although, 
the sign of the expansion coefficients stabilizes to 

all plus or all minr.s beyond a sufficiently high index m=M, 
the difference in the corresponding computed values of tKxjE) 
for tliese neighbouring values of E is found to be less than 
0(10“^^) for all 1x1 < x^^, r-rhere x^^ is the distance from the 
origin to the point in the nonclassical region at which 
t|^(xjE) begins to increase in magnitude. Saoisfying the 
virial theorem by the computed eigenfunctio'ns was used in 
the earlier literattire (e.g. Chan and Stelman 1963) to test fj 
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accuracy. However, ‘che fulfilment of the virial theorem is 
a necessrry lout cot sufficient requireraent (lOwdin 1959). 


I’he nOj?m of tlie confuted eigenfu.nct ion is. 


+00 


'tjL 

'2 ' f f (n) (n) { ^ _m+l _‘-2ax.‘~ 


r doc = I I j 


m 1 


'1 


2C 


doc , (2.19) 


"X 


'A 


where the rs.nge of integration is truncted a.t or^ — the point 

in the nonc.lassical region at which the comp^^ted eigenfunction 

\[jn(^^) tegirs to increase in magnitude. the value of if<^(x) is 

s^^ff iciently sma.ll lor |x| > x^ as dioctissed above and the 

conuribu-tion to the normalization from the rest o.f configura- 

bion space is estimated to be <0(10 ). '.Ohe integrals in 

eqii. (2,19) ftre obtained recursively starting from the 

1 

incomplete Gaussian integral J e dz (see Appendix B) . 

the plots of some normalized computed eigenfunctions are 
/ \ 

sho-57n in Pig. (II.l). 


I I • 7 Peature s of th e Method 

In this section some features of the method are seen 
in comparison with the other methods used for this eigenvalue 
problem. 

(i) It may be noted that no integration or diagonalization 
is necessary in this method which makes it attractive for 
the eigenvalue problems of the linear operators. 

* The vertical lines in the figure correspond to the classi- 
cal turning points. 
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(ii) In perturbotion theory uitli 
bills problem belonys to the singuJL 
■jie nriifoimi applicability of this 
its nonperturbcibive cha.racter. 


^ as the small parameter, 
rr" pert-urbation class, 
methoc" j?or all underline 


( iii) In Tiew of the becbion II. 1 an emp)ansion for the eigen- 
functions like (2.1) with a = constant may be cal.led as a 
‘fixed scale' exioonsion, buch am expansion is suitable only 
in a small re;ime of voilues of (njX) -vihere t.Pe scFile happens 
to bo close to the apimopriate ■value and it becomes unfavour- 
able in the other regimes of (n,X). 'Clie variational p.nd 
most of the numerical met]iods app3.ied earlier use the expansion 


(s:jX) 


- x^-/2 


00 

1 

m=0 


a X 
m 


m 


(2 .’ 20 ) 


in all 3oegime*s of (n,X). This expansion has a fixed scale 
a = ^ and is siiitable onlj" in the 'near hai’monic' regime. 

It is therefore not surprising that for hipheii o- X (nX» 
the eigenvalues could not be accuratel.y calculated in the 
above ¥orks. 


(iv) 'i'he need for introducing a scaled basis i-ras also 
realised by Keid (1970) who used the linear variation method 
for the pu-re q-qiartic osc'i-llator eigenvalue problem. However, 
in 9. variational framework the use of a scaled basis becomes 
intractably laborious for the following reasons. In a var- 
iational computation the first n (say) eigenvalues are 


od;h 
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obtained together. Since the appropriate scaling is dif- 
ferent for different n and X a single scaling is not suitable 
for the computation of alH n eigenTaluefj. A compromise 
scaling must therefore be used. But as n inc^reases this 
compromise scaling becomes imfavoiirable for more and mol’s 
eigenvalues. I'lie wo.yout is to compute ee.ch eigenvalue 
separately using an appropriacely scaled set of basis func- 
tions. this is intractably laborious in a variational scheme. 
In our method the use of an appropria'c;el 3 r scaled basis merely 
requires that a proper value of the scaling a obtained from 
fonnula (2.10) be used in recursions {2.3), (2,7) and (2.8). 
Since ea.ch eigenvalue is computed individually there are no 
carry over errors. 

(v) Computation ^.rith a larger basis is very simply done in 
this method bj'' continuing the recu.' sions (2.5), (2.7) and 
(2.8) for increasing m. In contrast a variational calculabion 
lutli o. larger basis requires integration and the subsequent 
diagonalization of a largo matrix which beyond a size is 
intractable, for insbe-ncej, the 10000th eigenvalue of the 
quartic anharmonic oscillator stabilizes to a 15 figure 
accuracy (in 3 minutes on IBIi 7044) at a point in the yecur- 
sions which corresponds to the use of nearly 17500 terms in 
the expansion (2.1). A variationa.l calculation of this size 
is inconceivable. 
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(vi) I’he method of infinite determinant for the eigenvalue 
pro Diems (Whittaker and Ij'atson 1927) wac used by IZerner (1951) 
and hisT-fas et al.. (1973). fhese e-ttempts ha,d a limited 
success. the primary reason for ttis is the use of a fixed 
scale expansion as discussed above. In the typical case of 
the quartic anharmonic oscillator (Biswa.s et al. 1975), only 
the lowest 8 oigenvEilues couD.d be obtained (using bho ezpan- 
Sion in the basis functions {x^ e~^ , until the numerical 

errors become too severe. Besides, in this '/ork, the pro- 
cedure used for evaluating the eigenv'i-lues consists of 
expanding the characteristic polynomials in powers of E and 
then finding its zeros. This is n-umericaldy inadvisable 
(see I'ox and Mayers 1968) because the uncertainties in the 
coefficients of the polynomials are hitphly correlated and 
the expanded polynomials with rounded coefficients are badly 
conditioned with respect to its zeros. Biswas et al. also 
carried out numerical investigations on the 'amount of 
nonnalization and the extent of orthogonalit3r ’ of the 
computed eigenfunctions in order to test their correctness. 

The overlap integrals required were evaluated by' integration 
over X from ~<o to +«o. ¥e have seen in Section II. 6 that 
an3r solution of the Schrbdinger equation *( ;E) -► +«> or - « 
as X ■►+«>, unless B is exactly cqua3. to an eigenvalue. 

Even a truncated expansion (truncated at a certain high 
index) gives a hiomp in the nonclassical region in the computed 


c 

1 

( 

1 

1 
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r 

i 

i 
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ei{^;enfun.ction, whore the actual ci^crif unction decays monoto— 
nicc'illy. 'Qie overlap inte£;ral£i evaluated hy integrating 
over infinite limit. s arc therefore ina.ccurate and the test 
of the eigenfunctions used hy Biswas ot a.l. is inconclusive. 



TA'RTn-; (ll.l) ; Str.biliza tion of the computed eigenval' 

(System’ .the Qu=.rtic . nharmonic Oscill- 
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I'ABLL (II. 2 ) z Comparison of our Results for B (1=1) 'witli the 

Results of Earlier Calculations?' 


Quantum ( X = 1) for the Qnartic Anharmonic Oscillator 

nuinber 


n 

This xTork 

Earlier works 



1 S' ' t l-( 



1-^387 1.87-78 (a) 

0 

1.392 351 641 53 c 29 

1.392 351 641 530 29 (b) 


« 

1. 392 350 653 579 1 ( c) 

10 

53.449 102 139 665 3 

53.448 ^ 0.1 6 (d) 

100 

1 035.544 183 138 91 

1 035.544 0/ (d) 

1000 

21 932.783 710 666 9 

21 932 '. 783 6 (d) 

1000 0 

471 103.777 790 8 O 9 

471 103.778 (d) 


( X = 1) for the Pure 

Quart! c 0sci3.1ator 


This work 

Earlier works 

0 

1.060 362 090 4G.1 18 

1.060 362 090 48 (e) 



50.256 25^1 516 7 (e) 

10 

' 50.256 254 516 632 9 




50.256 254 0 (d) 

100 

1 020.989 992 105 37 

1 020.989 99 (d) 

1000 

21 865.262 118 137 7 

21 865.262 1 (d) 

10000 

470 790.294 427 023 

470 790.293 (d) 


(a) Hioe and Montroll (1S75) - Eqn. (III. 9). 

(b) Biswas et al. (1975). 

( c) G-3i*affi et al. (1969). 

(d) Ilioe and Montroll (1975) - sophisticated hKB values. 

(e) Reid (1970). 
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TABLE (II. 3 ) i Eigenvalues of the Pure Quartic 0sci3-lator 

/ TT 2 4 > 

l,H = p + Xx j and the Quartic /mharDionic 
Oscillator (H = + Xx^') for X = 1. 

Quantura Pure Quartic Oscilla,tor Quartic Anharmonic Oscil- 
nuDiDer EigonveJ.ueo lator Eigenvalues 

ri 


0 

1. 

060 

362 

000 

484 

18 

1. 

592 

351 

641 

530 

29 

1 


799 

673 

029 

801 

40 

4. 

643 

812 

704 

212 

08 

2 

7. 

455 

697 

937 

986 

74 

8. 

655 

049 

957 

759 

31 

1-7 

J 

11. 

644 

745 

511 

378 

2 

13. 

156 

803 

898 

049 

9 

4 

16. 

261 

326 

018 

850 

2 

IS. 

057 

557 

436 

303 

3 

5 

21. 

238 

372 

918 

236 

0 

23 . 

297 

441 

451 

223 

2 

6 

26. 

528 

471 

183 

682 

5 

28. 

835 

338 

459 

504 

2 

7 

32 . 

098 

597 

710 

968 

3 

34. 

640 

848 

321 

111 

:> 

8 

37. 

925 

001 

027 

034 

0 

40 . 

690 

306 

082 

106 

4 

9 

43. 

981 

158 

097 

289 

7 

46 . 

965 

009 

505 

675 

5 

10 

50 . 

256 

254 

516 

682 

9 

53. 

449 

102 

139 

665 

3 

11 

56. 

734 

21A 

055 

173 

0 

60. 

129 

522 

959 

157 

8 

12 

63. 

403 

046 

986 

718 

9 

66. 

995 

030 

001 

247 

2 

13 

70 . 

252 

394 

628 

616 

6 

74. 

035 

874 

359 

102 

5 

14 

77. 

273 

200 

481 

984 

.0 

81. 

243 

505 

050 

767 

2 

15 

84. 

457 

466 

274 

942 

0 

88 . 

610 

348 

800 

799 

2 

16 

■ 91 . 

798 

066 

808 

991 

2 

96. 

129 

642 

045 

234 

1 

17 

99 . 

.288 

606 

660 

493 

3 

103'. 

.795 

300 

322 

273 


18 

106. 

923 

307 

381 

753 


111. 

601 

815 

045 

173 
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I’ablG (II. 3) (..Contd.) 


Qu-'-iittUD lure Quartic Oscillator Quartic AnliariDonic Oscil- 
iiumbcr Eigenvalues lator Eigenvalues 


n 











19 

114 = 

696 

917 

384 

585 

119. 

544 

170 

733 

050 

2-^0 

122. 

604 

639 

000 

999 

127. 

617 

777 

735 

355 

21 

130, 

642 

068 

748 

630 

135. 

818 

417 

325 

610 

22 

138. 

805 

147 

911 

395 

144 .. 

142 

195 

296 

398 

23 

147. 

090 

121 

257 

604 

152. 

585 

504 

205 

574 

24 

155. 

493 

502 

268 

682 

161. 

144 

990 

694 

513 

25 

164. 

012 

043 

622 

865 

169. 

817 

528 

001 

595 

26 

172. 

642 

711 

962 

845 

178. 

600 

192 

366 

876 

27 

181. 

382 

666 

185 

768 

187. 

490 

242 
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910 

221 
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199. 
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205. 

,582 

346 

604 

423 

30 

208. 

2 32 

539 

005 

144 

214. 

.773 

683 
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177 

31 

2-17. 

384 

261 
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103 

224, 

.074 
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052 

600 

32 

226. 
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4 81 

138 

233. 

. 466 

087 
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245. 
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050 

007 
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274. 
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.794 
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284. 

,068 
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581 

401 

291 

.723 

593 

051 

013 
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TABIi; (II. 3) (...Contd.) 


Qunntum number Pure Quartic Oscillator Quart ic Anliarmonic 


n 


Bi^envalueo 

0scil3.at 

;or Eigenvalues 

39 

293. 

948 

458 

266 

006 

301. 

736 

20 

351 

187 

40 

303. 

912 

066 

348 

384 

311. 

831 

516 

269 

701 

41 

313. 

.958 

030 

183 

978 

322. 

008 

069 

744 

845 

42 

324. 

085 

020 

992 

133 

332. 

264 

6U3 

530 

091 

43 

334. 

,291 

762 

334 

482 

342, 

593 

875 

832 

5.47 

44 

3^1. 

,577 

026 

891 

585 

353. 

012 

690 

233 

780 

45 

354. 

.939 

633 

506 

395 

363. 

501 

394 

863 

479 

46 

365. 

,370 

A 44 

467 

063 

374. 

066 

379 

800 

092 

47 

375. 

.892 

363 

004 

933 

384. 

,705 

074 

675 

721 

48 

386, 

.480 

330 

986 

517 

395. 

416 

946 

465 

263 

49 

397. 

.141 

326 

780 

674 

• :-06 . 

,200 

997 

442 

128 

50 

407. 

.874 

363 

284 

438 

417. 

056 

263 

284 

848 1 
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■l'ADLE (11.5) 


Eigenvalues of 

/ 2 2 
( n = p 4 - X + 

01 (n, X ) . 


tlie Quartic Anhfurnoiiic Oscillator 

/i 

\x‘') in Regimes oi Extreme Values 



0 

1.000 074 986 880 20 

36, 274- 458 133 736 8 


(near liana onic regime) 

(ne£.r qurrtic regime) 

louu 

2 154.242 0/15 2 52 21 

747 783.4-21 502 034 


(ooundary layer) 

(near auartic/VJICB regime) 




cmPTs±.. Ill 


Tia' G-EiTEiiAi hWiimmic osciiL/rion 


III 1 Introduction 

'i’h.e ci"env-'i,lue prolDlem of tiie /'onera.I anliarmonic oscil-- 
l;j,tor clesc^'ibed iJ3’’ the Hari'iltonian 

H = , (5.1) 


\rl,ioro X > 0, !i = 5»4, . .. and p = ~ i is considered in 

bln:, ch.ipter. Usinc the ,'calin, 2 , ar^inaents siiai.l.ar to rhose 
in (Jhnpter 1, it follo^'fs tTuat the eigenvalues of HCk^X ) 

=- + Xx^’’' are givcni "oj B^(k,X) = k ^ )y ‘''■fhere 

X' = . It ensures that the eigenvalue problem of the 

anliurmonic osciljator E(k., x) co,n bo corii')3.e te3.y described in 


box’mr of the I'educed Hamiltonian 11(1 s x) . 

A straifhtf oruard perturbative solution of this problem 
Lnns into d if fi cu.lt ies. ihe perturbation expansion for the 
eigenvalues in pouers of X is not convergent hut as3T33ptotic 
(bimon 1970). Ibo coefficients in the perturbation series 
grow vG.r3- fnst and the construction of v. .rious Pad6 approxima- 
n.ts bo COL) e extremely involved, the numerical resul'co of 
Craffi et a.i , (1971) for the octic a}Tharmonic oscillator 
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( p,=4 ill- (3-1)) suggest tliut Pads approximants do not converge 
to the exact eigenvalue. Ilouever, mixed Sorel-Pade method 
■(■’■ere utilized to obtain a few eigenvalues. iTon-perturbative 
calculations have been relatively more successful. Biswas 
et al. (1973) extended tl' j 'Hill dDteriiiinant ' method, used 
for blie qaartic a.nharmonic oocj.lla,to-.- problem, for this case. 

p 

m — 

'L'liey ii.'.:ed the basis functions {x e ' ioi" the expansion of 
the eigenvalues for the sextic and the octic anharmonic 


oscillabors ( p,=3 and 4 respectively in (5-1)) for values 
ci X i'n. the range 0 < X < 100 for n=0 and u < X < 10 for n=2. 
However, the accuracy of the computed eigenvalues in their 
work reduces significantlA?- as one goes from the qus,rtic to 
the sextic or the octic anharmonic oscillator, besides, the 
evaluabiori of the eigenvalues gets confined to smaller regime 
of (n,x) values. Lakshmanan and Prabhakaran (1973) obtained 
comiclassically an asymptotic expression for E^^‘'^^(x) in the 

p,=3 case, truong (1975) used Weyl-quantization prescription 
to study the nextio anharmonic oscillator eigenvalues bxxt no 
new re.sults were obtained- For sufficiently large n, the 
111 KB apiorozima-t ion method has been used to obta.in approximate 
oiigenvalues . 

■fhe most recent and extensive work on the genera.! anharmonic 
oscillator eigenvalue problem i: due to hioe, MacMillen 
and Ffontroll (1976). fhey distinguished two limiting regimes ; 


i 

i 


I 
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oi valAieo of (n^X) for the eigenva.lues enalogous 

to the fiimi.Lar diotinction made for the Cjuartic ahharmonic 
ofjcilliator oigeavaliiGo. In one regime the ener^^^y eigenvalues 
diffex slight l_y from the harmonic oscillator leveJ-s (the 
hoar bo, rmonic' regime) ; in the other they differ s3.ightly 
from the pure 2j.'.--ic oscillator eigenvalues (the 'near pure 
anhairaonic' regime) . the above tvo limiting regimes are 
,separa. bed by a regime ca,lled the 'boundary layer' in uhich 
che energy eigenvalues are not 'neardr- ha.rmonic' or 'nearly 
pure enho.rnionic ‘ 0 Hioe et al, used different formulations 
of the eigenvalue problem in various regimes and constructed 
several simple formulae uith different ranges of validity. 
Using hargi.iann representation, they developed numerical 
a3.gorithras from which tho energy^ oigenvad-ues in the smald. n 
regime may ho computed. The algoirithjn is simi3.ar to that 
fiovo3 opod foi." the quartie anharmonic oscild.ator eigenvalue 

px"Ob3.em . First few Gig;onv3.1i''es xrere thus computed to 5“6 
significant figures foj- t]K-5 sextic and the octic anharmonic 
oscil'Lr bors, loi’ variowis Vwilue of X. Loxfever, tho size of 
the determinants, required iox the computation of the eigen- 
vr.,li>os, increases rapidly eith p o,. n, making the evaluation 
of h.ighsr cigenvad.ues laborious. 

I'Je sho-.’f in this chapter that the method described in 
Chapter II may be extended to solve the general anharmonic 
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03 ci.l lat 01 ” eigenvalue problem and eigenvalues of arbitra“ 
rily high accuracy are obtainable in this case in al3- regimes 
of (n,x) as for the rj_u.artic anharmonic osci.llator problem. 


1 1 1 . 2 f]ie__ Sca ling Por mula 

iiie eigenvalues and the eigenfunctions 

genera3. ahnai-’monic oscillator (H = + 

T?- + are the solutions of the Schrbdinger equations 


[„ 1™ . + x2 Xx^^] Ip ) , 

(3.2) 


CIX 


n 


T^'itli the boundary condition (xj x) -*■ 0 as x-^+<». he 

write eigenfunctions in the form 





00 


I 

m=0 


a 


m 




(3.3) 


where c; in tlie scaling constant, A formula for determining 
the appropris te .scaling a for any a i,n obtained along the 
Scirae l.ines a,o for the nuartic anharmonic oscillator (section 
IT,1). I’he exact n-th eigenfunction has n zei-os in the ' 
region of osci.llation which, for sufficiently large n, is 
v(E 2ii)y j^j.l/2p.^ yjQj approximation 

5 i ( 2 p)(^ ) = C (n + i) ' 

XI ^ 

. -nn .r 2TlltlT /^,lJ?+l)’ 

Hence, the region of osci.llation '^X m / (n+^J 

Eor an effective expansion bhe region of oscil-lation must 


CEh,^ 






^ 8469 . 
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include tlie span of lirst n (or a numlDer- proportionaJ. to n) 
lasis functions. The appropriate scaling formula is thus 
obbaincd Top setting the region of oscillation of the nth 
eigenfunction equal to the width of the nth basis function 
(''' 0 :"'"^'^^ , which yields 

u-1 1 

■v (n + i) ai xii+i . 


j.ito 1 -above scaling formula is not expected to be good when 
n nj\d/or x i-' email, in Tnew of the yKB approxiriiation for 

iil used in t]ie derivation. However, for smsi-11 n, 

II 

f;r,ifil.l X the scaling must approjch 1/2 — the appropriate 
scaling; for the harmonic oscillator. Hence, the scaling 

lormula valid in al'i regimes of (n,x) and for any p is 

lirA 

a(n,x) = w i (n + i)^'*'‘'*‘X (3.4) 


Thcj following regimes 'may be distinguished according to the 
o.bovu scaling formula; 


(n + 


ihhl .1.. 


1 

2 


(n + 


Ikrl 

1 ^ p+1 


1, 

p+l 


>> 


1 

2 


(n + 


ikzi • 

Ij p+1 ~ ■— 


is the near harmonic regime. 


is the pure anharmonic regime, 


is the boundary layer between 


the above two regimes 
The value of the combination (n + 2 )*" ”X determines 
the regime to which the eigenvalue 


belongs. It may 
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be seen from above that the range of (njX) va3.ues for the 
’neo.r harmonic regime* dimnishes on increasing \x„ Hence 
the eigenvalues of general anharmonic oscillators, obtained 
from the methods which use a finite term expansion for the 
eigenfunctions in the harmonic oscillator basis { ^ , 

deteriorate in accuracy on going to oscii.lators of higher p. 
fhe use oi basis, appropriate 13'- scaled according to the 
regime, loads to a uniform treatment of the anharmonic 
oscill.itors eigenvalue problem for all u and in all regime 
of (n,X). 

Ill . 3 Method 

The expansion (3.3) on substitution into the Schrbdinger 
equation for the general anharmonic oscillatox- (3.2) yields 
the following (p+2)“term linear rec^■lrrence relation among 
the expansion coefficients | » 

(m+l)(m-i-2) a^_^2 (E~4cm-2a) a^ + ^ ^ ^m-2p ^ 

(3.5) 

The even and odd parity so3.utions are obtained respectively 
by assigning the initial conditions (i) aQ=l, a'-j_=0 
(ii) a =0, a-|=l. lie divide the recursion (3.5) "by (m+l)(m+2) 
and rewrite it in the folloi;fd.ng notations 

^■m+2 ^m,m % ^m,m-2 ®'m-2 ^m,m-2p ®'m-2p 


= 0 , 



vrhere 
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(E-4am“2a) 

J- 

(m+iTTm+oT ’ 

^m,m-2 = 


T^i-iTr35-27 ^ 


(3.6) 


“ Tm+]7Tm+2T ' 

For BGlf coniiistoncy the determinant 


A(E): 


^00 

dao ^22 


0 

1 


^ '^44 


U 

1 0 


0 d o-Od „d 10 

m,m“2|j, m,m-2 m,m 


=0. (3.7) 


The eigenvalues E„^^*^^(x) are the roots of this tra,nscend 3 ntal 
equation. To obtain the roots numerica.lls'’, I’o denote the 
detorminont formed by onaittin;-'. all ro^^s and columns beyond 
d^ ^ in A (E) as determinant 4^+2*^^^ 

expanded into doborranants of lover orders. It 3 i^ieldD the fol- 
lowing ( n+2)-term rocuin'encG relatiOj.: among lAjm-. 




d (h) A (E) + d _ A ^(1) 
m m>m“2 m“2' 


+ (-l) 


p-1 


d . A ^ (E) = 0 

n,m“2p u-2|i^ 


(3.8) 


Hence the values of the determinants upto any order 
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’Tiay bo determined successivel^r in terms of Aq (or A-j ) ‘^'ritli 
tliG help 01 tlio ci.oove recursiono flie corresponding^’ zeros 
of > ‘ ®’^'^'bilize to faa eigenvalue for large 

m pi’Ovided the scalin{;' oc is appropriate- The procedure 
for the actual numerical evaluation of the eigenvalues is 
the !-,.ame an ir, described in section (II. 2)- The PTouton’s 

iUOGiiod,, which in used lor obtaininr the roots of A^ (E)=Oj 

in 

roqniros the value of the derivative A' (E) also. The deriva- 

m 

tives Aj^(E) can also bo evaluated recursively with the 
ho.l p of the recursion s 




in,ni-2 ^m-2 


a: ,(b) 


C-O'"-" 4h2^(E) = 0, 


( 3 . 9 ) 


cbt ined by differentiating (3*8) with respect to E. The 
I’OCTirsions (3.8) and (3.9) arc continued on computer until 
tlio corrections for the required root of {A^^(S)} for suffi- 
ciently l.'^rgo m stabilize to a proscribed extent. The 
eigenvalues accurate to 15 significant figures may thus be 
compubed in on3-y -'i or 5 itcrctions starting from rather crude 
initial ostinatec. Tho initial estimates for oigonvalues for 
largo n (and not too low x) are obtained from the corres- 
ponding; VIEB approximation formulae and fox' lo”^ n tiiejr aro 
obtained by ovalurting a sufficiently large order determinant 

A (1) from recursion (3.8) at various E points. 
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For obta-inin.'i: eigonfimctions, 2 ic sec eatial to tlie 
comprtod Gigonve.lues in f'.u recursion (5.5) and the expansion 
coofxicients arc evaluated rniccessi vely 'i'ho 

roGultin^' coefficients {a^^(E)} provide a convenient repre- 
sonbci-tion for t]io corrospon.din/ eigenfunction tliroou' li the 
cxpo-iir ion (5.5)* -fac asymptotic bohavioeir of the solution 
of the Sclir'ddin/ or ociua,tion (5.2) v exp (j-Jx} ^"^"’‘/[.i+l) , which 
.requires 


a 








m-*- TO 


( 5 . 10 ) 


for Iho coefficients in bhe scries solution (5.3). It is 
•.■ 0,1 during the computation of {a^(E)} that the above 
asymptotic dopondence is satisfied actually;, rdiich ensures 
u cie creasing (in magnitude) set of coefficients for suffi“ 
cjcjitly -lai'go m. It is ao'w possiblv. to obtain accurate 
cjgonvaluos and oigonf unctions for the general anharmonic 
oscillators for ary value (n,X). 

1 1 1 . ^ 

ho hiive obtained the eigenvalues of the soxtic 

oof 2*28 

„ p2 + Xx°, ^ >0) and the octic (lI - p •+■ x + Xx , 

X > O) anharm.onic oscillators for various values of n 
and X* eigenvalues are computed using values of a in 

the range: 
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a( XI, X) 


= I + (1.5 to 2.0) (n + , (for the 


sexibic anharmonic oscillator) 


(3.10a) 


a(rL,x) = ^ + (2.5 to 3.0)(n + y) , (for 


'cne 


octic anliarmonic osci.l3-ator) . 


(3.10b) 


Olio ran^'C of a above works for any value of (n,X). 

I'lio oiyuiivaJ uoc , thuo computods are prosoiited in Cables 
(lll.l) .'.rd (III. 2) for n = Oj, .Ip 2p .. 10, 100 and 1000 
and j.or different va.luos of X in the remge .00001 <X< 40000. 
Choy arc cvaluatcKl to 16 lif^Tor-os and blion arc rounded off 
bo 15 figures fo-o tho '.Cables, fhe 'labJ.os (III..1) and (III. 2) 
inciLudo vodtios of E^^^'“’'^(x) for the (high np lou x ) region 
o;. the boundi'.r 3 r lover. Ihc compubation of the eigenvalues 
j.n this region ic found to bo tb.e most difficult in the 


ova lior literature, Itirthcr. the va.lues of E 


ji (2 |J,) 


n 


(x) for 


bhc soxtic and bho octic anliarmonic osci.llators to this 
■ iccnr. cy jn .il" rcgiiaos of ( i, X) are ropoited here for the 
iir,ob time uxio'.,ing' most oxtenrivc tables for E^^^^^(X); 

|.i = 3, 4p aro due to Hioc et a'’. (1576) who have presented 
tho V. '.lues or first six eig,cnva.lucs to G significant figures 
for tho ooxtic anharmonic oscillator and of first four eigen- 
va.luos to 5 oi{;nificant figures for the octic anharmonic 
oficillator. Iho possibi.lity of treating all higher order 
gonural anharmonic oscillators oigenva.lue prob3-cm by the same 
toebniouo for any va.luc of (iipX) is thus established in this 
chapter. 
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CEAPIBE lY 


TI-I3 DOUBLE MIDIKUI'f uSCILl-ATOB 


17.1 r p due t i on 

¥o consider in Ihir. chapter the ei,penTolue problem oj 
the double minimum oscillator (d.ir.o.) desc.ibed by rhe 
jiamiltonian 


n(i,x) = p^' - . X > 0 . 


(4.1) 


iho potential function of a d.m.o. has tuo S 3 "mmetric poten- 
bL.J. wells separated by a bar-rier. A feature of its eiaen- 
voluo problem is tbe bunching of the lower eigenvalues in 
pairs for suff icien'cl 3 r largo separation between the two wells, 
iho d.rn.o. models some interestin,': ph 7 /i ical problems. 'Ihe 
vib ’.tional specti’a of some molecixl-as possess t^ro pax’allel 
tjrpe nearly superimposed bands, a phenoDienon which mo.y be 
directly/ related to the eigenvalue spoctruii of the d.m.o. 

I'he commonly knovm example in this regard is bbe inversion 
spectra of the ammonia molecule (see Doin,Lison and Ulolenbech 
1952). Besides, the potential functions of several hydrogen 
bonded solids are found to possess two minima in the region 
available for pro conic movement (see Synder and Ibers 1962, 
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k->o.'iOj?.i.-ii and Horiii{j, 1362) . llie most recent appearance of 
tiTx.'j nooel is in the spontaneotic; S3n?netr/ breaking (Polyakov 
1977) . 

Bie coiTiputa bion of the sp3.ittiiif' betveen the ei;j;envalues 
i-Ormin^' pairs lias ijccn the subject of cor sic’crable interest, 
I'hc '■^plittiiig depends in gonoral on the separation betiiveen 
blie tuo wells and the nature of harrier betxireen tliem^ The 
d.m.Oc with pobential function V(x) = ■^■'k (Ixl- a)^ is exactly 
wolvo-ble (Merzbacher 1961) o-nc' an expression for the splitting 
botwcGii the tiTO lowest oiijenvolnes is Af v (^ka^), 

honnison and Uhlonbock (3.932) obto^ined the splitting in the 
'.KB ax)proximation and then compared the hhB va3_ues with the 
OX: ot vo,.lue:'' for a d.m. o with potent ia,! function fomted by 
joining two equal parabolo.o with a straight lin-^. The I'EB 
v.'tlues for the splitting arc found fair3.'y accui’afce for large 
sej)P.ra fcions botwcun blio t'/o pax’r.bolae. It is inberosting to 
note that bhc '/iJi upiiroximr/bion is applied here for low n. 
j.'tio usefinlnof-.s of t]''iOHO results is, liowever, limited to some 
extent, as bhe potontiad. fmebions used in these works are 
non-aTialytic . Harmony (1973.) trc-sbod tho chm.o. problem 
Vi i a hs.rruonic osciHo.tor approximation and obtained zero- 
order and first-order expressions for the splittings 

■L'he pL-rturbation expansion of tho energy eigenvalues 
of tiic d m.o. described by (-t.l) in power series of 1 is 
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uoii-'CC^nvorgeiit (Simon 1970). Somorjc'.i an^'’ Hornir, (1962) 
oOtaiiLocl uuiacric-‘'Llly o, oiiercy ei£!,GnvalviG& ( bo 9-5 

sj.. .nil icai'ii; fij-'nireo) for the d.m.o. iritli Ilarfiiltonian 
][(kjX) - p'" ■■ -1- Xx^ for five different paiig of values 

oi' (ksX)- fko c -Iculations were done Ly expanding tlie 
oi/ on j.anc cions in tlic ho.imonic oscilla,tor Laois functions 


and diayona.liziix^, tlio secular determinant formed, le oLtainj 
in this cTaaptor^ accurate eigenvari-ucs and eigenfunctions of 
(4.3.) using the rnotlod described in Cha-ptor II. A VO 
expression for the splittin/. is e.lso obtained for this 
X^robJ-em and. the VKB vri3.uos are compared with the cones- 
ponding accunsite vs lues for the split bin® dor various values 
oI X. 


I . 2 .i|i ifeji v. t Itip s_ 


,:)Cli3''ddinger equition for the d.m.o. (4.1) is 




2 


dx 


i Ix^] 


E^(X) ^JxjX) . (4.2) 


fb>; ei.pcnf unctions are expanded as 


,{,(xix) = I ^ 


m 


(4.9) 


m=0 


vrhich on subsjtiitution into (4.2) yields the following 4 -term 
recurronco relation ^ 
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ni+2 nijin m in,in-2 jd-2 ° 


m;ra“'4 ‘^n“4 


= 0, (4.4) 


where 


d. 


m 


- .Ill,’; 4 ct-1“2g;) , 

sTP. ^+1)(ei+ 2') ’ ''rniii-: 


4a + 1 
2 = ’ 


^ro,Ei^4 " rn^^r)(u+2T “ 

■fho motliod of corjputino; tlio oigcnvo.luoo from a rocTjrsion of 
'Glic type (^ 0 ^) lr.io boon d escribed in Chapter II. fhereforOj 
bb- oi/'envc'-lues of the d.iiuo. arc obtainable Tritli uniform 
Mxou,:o,cy in all regimes of (n,X). fhc charQ,ctcristic 
i)uuc}iin; of tlio eigenvalues in p)' 2 -irs occurs for smo.13. n 
‘■'''If' I or sufficiently large separation bot^-’'eon the tf’jo wellB. 
I'JiJico blie sepe,ration butweon the tfro ifells the 

roglon (;f intoro.st for the present problem is the (low n, 

■lof- X) rogimo, lor (high n, high X) regime the eigenvalues 
are iiee.r pure quartic. Ifo have therefore compirbed the 
cige/iv/i.lucs in bho (low n, low X ) i-ogimo and. the results 
ai'e presented in fable (IV. 1) ior bhe firsb eip,ht oigjon*- 
values for viluos of X in bhe range C'.0.1< X < 0.20. fhe 
values of the scaling a used in those computations lie 
bi.'twecn 0.5 and 1.0. fhe eigenvalues prooented in the 
.'.‘able (IV.l) arc with respect to the bottom of the potential 
wo3.1.y ab zero cncipgy and is related to l^^Cx) by 


Enb). = I' + 


(4.5) 
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■'‘'rtiere l/4X is the depth of the potentia], rell. The numbers 


e^( X) are positive definite and provide a direct look at 


the variation of the eiAenvalues ’.rith x 


JJ O-i 


0 .. the 


' iei-)r'.T€it±on between the two vrells increases and the prohahi“ 
lity of penetration throu^’h che harrisj"' approaches zero. The 


e igenv. .lue s e 2 


and e 


2n-rl’ 


therefore,, becoiiiO nearly do^ably 


doyenorato for sm^ill n. For instance , £ 2 ^ ^2n-i-l 


(n = 0,1;, 2) are found close to each other to at least 14- 
si.^ nificant figures for X =0.01. I'lie expansion of the poten- 
bial function of the d.m.o, about bhe minima of the well 

for X ->■ 0} therefore the lo^rer ei,-" envalues a,ir2(2n+l). 
'I'he numerical results confirm this observation. 


The eigenvalues of Il(k,x) = p^-kx^+Xx'^' are obtainable 
from the eigenvalues of 11(1, x’) using the sco.ling relation 

Eglc.x) = 

where X' = k • 


IV . 0 Th p, ]J ig^. kprpulp, for S pli tting 

iho splitting; between the pairs of loxrer eigenvalues of 
a syminetric cl.m-o. in The hKB appro ima’cion is ^,iven oy 
(Landau and lifshitz 1965) 

,Xq 

^J,KB ^ a. exp r- ! IPldx], (4.6) 

where = h- P~^ and ±Xo, are bhe four turning 

Xq 



73 


points. ITlie derivation of tbe above formula assumos small 
probo-bility of penetration tliro^igb the barrier. For the 
cl. 111 . o, deocribod b^/ (/!.3.) tbe tuT>ninr points are given by 


— q”* (7“"'^) 5 ( l+n) 


(4.7) 


“ ^n 4x 


-^n mean energy of 


•0 “ 2a 

nhere n = '(iAxe^) 
tbo ttfo eigenvaluer; forming a pair. The inte^.rals involved 
in (4-6) may be expressed in terms of the completo elliptic 
intograls K(k) and E(k) of the first and tlie second kinds 
respectively (G-radsbteyn a.nd kyzhik 1965) 


+ Xq 

[ [Pjdx 
■^0 


J 


~\r2 


[{x2 


0 


2w 2 

X ) ( 


ax 




and 


2^1/2 ^1[(x 2+2,2^ -(t) - (x^-x^) K(t)] (4.8a) 


"1 -1 n 1 fl 2 2 n/ 2 .2xn~l/2 

.-L ^ -J- r r r ^r ^ dx 


./ = r [< 

^"0 0 


, 1-Ux 1/2 

miere x = (-^^ ) . Q. 


- K(q) 


(4.&b) 


.2u xl/2 


i.'hus s, 


AS 


USB 


2.i(.l(i+:^J., 

i:(q3 


1/2 


exp 


til, 

5 X 


( H-u) 


1/2 


mt) 


Xilx( t)) ] 

(4.9) 


fho values oi the splitting AE are calculated from the above 
formula, for variou'^i values of x for the lowest two eigenvalues 
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and arc compa.red i-rith tlie corresponding accurate values in 
iable (IV, 2) „ lie mean energy Ii° required in (4.9) is 
evaluated from iallo (Iv.l), The UiCB values for splitting 
aro surprisingly good, 

Poj’ Gii'cO.l X a simple anal'/txo e.piproniiiate expression 
for XE mu-y iiov irj ohtainod xron (4.9) using the following 
oxp. i, Visions for the elliptic integrol.s ( Cri’adshteyn and nyzhik 
1965) : 

L-(k) -= I (1 -1- -4" Ic^) + O(lc^), 

Edc) = I (1 - ™ k^) + O(k^), k 0 (4.10) 

and 

K(k) = In + I (In |t - 1) + 0(k'^ In kO, 

E(k) = 1 + I (In "T “ f) + G(k'^ In k’). 

k' = /( l“k^') , k Is k' -»■ 0 . ( 4 . 11 ) 
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Fig. IV-1 Double Minimum Oscillator Eigenfunctions 
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.LABjjji, (IY,2) Comparison of tiie WJvb Tallies fox* fhe splitting 

with the corresponding acciii’ate Talues 
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CHAPTER V 


THE TRANSITION MOMENTS 

Tlie tranrjition loomonts betxreen the anlir-rnoni c oscilla- 
tor onej. .".V Gi,.'_eTiGt.i. ces j ii> and ln*> are blie raaci'ix elements 
< n jn ‘ > . The best knoACi estimates ( ClMan and Stelman 1963i< 
.Roid 1970) oi the branoition moments were obtained for the 
pure ijuartic oscillatoi’ nsin.^ variational ei^''enf'unctions , 

Tig vari. .tional oipe/ilunct: ons are Icnonn to be much less 
accurate bban tlie corresponding oigenvaluco and are unsni’caole 
■.’or the computation of the transition moments of high accur 8 C 3 r. 
T.o .'cciravte evaluation of the eigenvalues and the eigen- 
lunctionrs in this work makes it possible to obtain accurabe 
trunc.icion moments. Eurther, the matrix elements <nlx^ln'> 
satisfy an exact linear recurrence relation in tie in(3ex k 
''janorjoo 1977)" The recurrence I'el.vbion makes possible the 
cjvnlu-’ tion of all hipher moments in terms 01 lot'cr nomenxs and 
oi. ■ onvalues j, x>ritl'.ou'c inbegivioion. 


V . ]. Th e i' ® cuiT Giicc felat i__9n__ f 

* 0 consider bhe class of llamil conic;..nfc. n — p -r Y(x) 
where V(x) is a real polynomial function and Kln> = E^ln>. 

Then 

<nl[H.h]!n'> = 


(B^ " <nlhln'> 


(5.1) 
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k - k“l 

fOj' rin arbitrary 0}^)erator V. On settin^c; = r. and r" 'p 
r.nicceGGively ir (5.1) and usin;, the &c>Lr'd dinger equation, 
f ol.l o .'ing re Ln.t.T one are obtained; 


(2n “ En 


,) <njr '|n'>= • ■k( k-'l) <n jx^'" ^jn’> “ 2ik<n }x" pjn' >, 

(5.2) 


^^n ■' ^'n'^ <nlx^'“^pln'> = -( ’.O'l) ( k-2) <njx''‘' ^p|n'> 


2i(l:-a) <nlx‘"“^ln'> 


+ 2i(k“l) <n Ix^^‘"^V(x)| n’> + i<n }x"'* V'Cx)! n' >, 
whore V’(x) = 3V(x)/3x. ( 5 , 5 )) 

III imino, ting the nv^trix eleraents of x^”^p and x" -^p from the 
above two relations, we obtain 


41c(k-l) <n|x^“^Vln*> + 2k<njx^“'^ T’j n'> 

= ( I ■; -lil ,)^ <njx^ln'> + 2k(k-l)(3 
'' n n ' ' ' 

-h k(k-l)(k-2)(k-3) <n|x"^“'’-in'> 


-;-E ,) <nlx^" ^jn' 
n n ' ' ' 


(5-4) 


> 


For a poljmomial pobential 
to the required reciirrence 
matrix elements required t 
polynomial potential V(x) 
involves the same elements 


Vlx), the equation (5-4) reduces - 
relation, flie number of initial i 

o start the recursion is la for a ^ 

I 

of degree 2p.. Iho recursion | 

of all the matrices. fhus the nn'-tlj 
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cleiiiont Oji the loweet \i matrices and the tvo eigenvalues 
v'.no E_ 0 arc suilicient to determine l;he nn^-th el'jment of 


n 


n’ 


X for any k* It may be noted that '’■fthcut tie above 
rucLirrJion one noecls cJ.l the matrix elenients of the matrix 
X be dobermine a element of In tlie coee of the 


cuartic anliarinojiic (potential funebion Y{x) 




+ xx'^) 


rind tl'ie pure quartic oscillabor (V(s;) = X:t ) the rela-cion 
(‘' 1 . 4 ) yields tto follovfinr necmreiicc relabions re£?pectivel 7 


4'..(lr-rl) X <n|yh2ln'> = <nlf''|n-> + 


2k(k-l) ) < n > 

+ k(k-l)(k-2)(k-3) <nlx^“^ln’> ^ (5.5) 

id 

4]'.(k+l) X < n n' > = <nlx^^ln*> 


+ 2k(k-l)(B^+E^,) <nlx^“2ln'> 

+ k(k-l) (k-2)(k-5) <n|x^"'^ln'> . (5„6) 


Thun all higher moincnts for- any particular transition of 
the cuartic anharrionic 03.‘ the pure quart ic oscillator may 
be obtained in terms of tlio lowest non-zero moment of that 
transition, fho initial requirements in these cases is 
effectively reduced to one moment because of the even 
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fj 3 ni.mo‘try of fio Ha-miltoiaiarij, the other moment bein-^ zero^ 

In I'able (V.l) tlio ra.tio of the Ul element for the matrices 

jc, ^ fo’.’ the nuartic anbari.ionic and the pure 

quii^'tic oncili- '.t oiT nr‘o tabnlatod. fov lauj^p k the reciir- 
ni ons ( 5 - 5 ) anc! ( D . 6 ) yield 


X <nlx^''ln'> /<nj n'> '-vk" 


k -> <» 


(5.7) 


which doscriboB Ihc asymptotic behc.-viour of the raoments in 
tl'UBO cases. 


v . 2 C omiru t atp.p;i p_f the, iransit ion Momen;;^, 

The lovrur moments between various pairs of tnc 
anhai-monic oscillator eigenstates may now be computed in 
Ihe following manner. The expansion coefficients {a^^Cl)} 
for the required eigenfunctions are evaluated by substitut- 
ing’ the corresponding computed eigenvalues in the recursion 
lor {a (E)}. h’'e include as many number of coefxicients 
in the expansion of an eigenfunction as were required for 
obtaining the eigenvalue stable to 16 significant figures. 
The range of inbogration for the evaluation of the transi- 
tion moments is trmcated at x = the point in the non- 

clasr ical region at which the computed eigenfunction just 
begins to increase in magnitude. bJince the eigenvalues 
used are accurate to 15 significant figures, the computed 
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cii^'enlunctiono roach estremely small values in the non- 
oT ifesical region before they finally start incroai,sing in 
laagnitude for large | s| (see section II.'o). The contribu- 
tion to tbo transition moments from the rest of the 
coni i;''ura,tion spoce is estimated to be <0(10'’^''^). 'fliuS;, 




<nlAln'> = z Z ai") a / ' e-{f 


3 3 


iiaoro |n> = Z a 

j 


“X 


■A 


(5.8) 




1‘ r4.iic^ N 3.rc *fcli6 13 ormaliza'tion cons'fcs.n'ts for tlie rospsefiv© srfcafe 

l‘\o intu{j,ralG involved xn (5*8) may be expressed in terms 

1 2 

of the integrate ^ dx, (s = integer), 

which satisfy the recurrence relation (see Appendix B) , 

(2S-1) l2c,..2(P) " (2s+l+2f3) l2s(P) + 2p I28+2^P^ = * 

(5.9) 

The actual proceciure for the evaluation of described ’ 

in Appendix B. The computation of moments therefore i^e quires [ 

i 

no integration. i 

2 I 

The non-zero matrix elements of x and x in the lowest j 

[ 

ten eigenstates of the quartic anharmonic and the pure quart ic | 

oscild.ators for X = 1 wore thus computed and are presented ■ 

i 

in Tables (V.2) and (¥.5). Each moment given in the Tables f 
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(Y.2) and (¥.3) lias lioon clieclrod by varying; a in the 
ap].iropriato ro.np;© and is cl.aimed to 'jO acciinatc to sJ-1 12 
-Cir-iiren yivcn in tallies, The transition moments for the 
quortj c anhnrmonic oscillator o,re roportod for the first 
bimo in this nook. Amoiiq bhe oirlior libera taro only a 
fen nor!“'ir.oro r; onicnts for fio pure ouartic or.ci] la,t or n'ei'‘o 
evaluated to some accuro-cy (Chan and Stelman 1963? fold 
1970). ihc coiv’os pond inf prosciib values for the pure 
quartic or-ici]..lo,tor moments arc ab .lease 5"6 figures more 
accurate, further, the recursions (9.5) And (5.6) give 
a] .1 the hif her moments to bhc same accurs.c 3 r ao of the 
3.ui'ust non--zcro moment for tht.rb traiisiti03a, witnout inbogra- 
t i on . 




Tho work presented in this thesis forms a part o± a paper 
entitled 'Tho Anharmonic Oscillator' accopbed for publica- 
bion in tho Proceedings ox tne Poj'al Soi.,iety. 
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i'A}3IiIii (V.l) ; The ratio [< 1 |x ' |0 >/ <l| 0 >] for the Quartic 

Anliarm or i c and the Pirre Quartic Oscillators 
from recursions (5.5) an''"’ (5.o)« 

__ [<1 |x"^l0^y<l|xj 0>] 

The Pure Quartic Oscillator The Ouartic Anliarinonic 
(X = 1) Oscil l ator (X= l ) 


1 1.0 1.0 



0.957 97 c 052 702 071 

V...825 567 331 595 326 

13 

1.361 642 584 655 12 

1.073 5 u 6 435 204 95 

7 

2.604 OGl 403 536 68 

1.862 740 950 928 00 

0 

G.086 204 025 521 49 

3.906 125 267 322 92 

11 

16.628 157 346 910 1 

10.048 610 500 904 6 

15 

51.631 352 902 97^/ 5 

20.965 390 251 421 1 

15 

.170.629 0 O 6 G.11 098 

93.500 722 956 452 3 

17 

678.615 hi 2 628 609 

332.840 550 182 597 

15 

2 799.464 881 619 44 

1 291.293 351 789 83 

21 

12 /!30.u58 107 522 6 

5 409.354 773 052 46 

2 5 

50 981.717 084 750 6 

24 284.238 551 754 3 

25 

297 337.097 04 I 684 

316 109.407 781 311 
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APPLiroiX A 


describe liere tLe procedure adop'ced for obtaining 
initial ectimates of tbe eigenvalues in the '‘bouncls.rp layer' 
WuxCr n is large. Let us sr.ppose that ube desired eigenvalue 
is (X = Xt^) > uLere (n, X -j^) Hee in the 'boundary lc,,yer*. 
iin initial estinate for ib is obtained as follows; 

( i) •■:e first obtain accur:-te eige3nv3j.ue (x = Xq) > wliere 
(n, Xq) lies in t]ie piire o-nharmonic region and the initial 
ostiina.te for it is obtainable from the ¥KP formula. 

(ii) file desired value X =Xi. is reached through a seouence 

iJ 

of intermediate values { X } s i = 0^ 1> 2 , ... . The initial 
estimate for (x obtained using accurate values 

of (X = Xj_) snd the I'aylor series expansion. It gives 


E, 


initial 


3l 


n 


*'\+l^ ^ ^i^ W-rl"’ ^i^ a 


n 


ax 


X . 

1 


i=l,2. 


vheue 


aE 

ax* 


n 


X = X . 

1 


^ n ^^ i-1^ 
Xi “ ^i-1 


i’ho values {Xj_} are chosen sufficiently close to each other 
so that E (x • , n) ciay be computed avoiding .-jumps to 


E 


h+1 ^^i+1^* 


( iii) iy refined to 15--fignre accurate eigen- 


va.l 


Lue B (X. .) by the method used in 'bhis thesis (Section II. 2) 


(iv) Steps (ii) and (iii) are continued till the value X=X 


b 


is reached. 



APPEHDIZ B 


I'hG I’eciixsivo evaluation of the integralo defined, by 

l2s(M 




(B.l) 


U 


i, conrjidereci in tliis Append On integrating, (Ihl) by 
parbKj, tre obtain a fol.Lovinp inhomo i^er-eoaio recurrence relation 


2p l2.s+2^^^ " 


(13.2) 


‘.evrl bin/]: (B.2) on replacing; the index 'o' by 's-2 


?P l2s(p) ~ (2G‘'1) l2s-2^P^ 

find eliminating, the inhomo{-:eneouB part from (B.2) and (B.3)» 
one obtains a 'i'--tern homogeneous recu'crence relation for 
123(13) : 


■2s+2' 


(p) - (2s+l+2p) ( 2S“lj l2s-2^^^ 


(B.4) 


fhe inbegralB l 2 p,(P) ™sy therefore be computed for any va3.ue 
of s by siiccesiiive app.licatioii of the relation (B.4) • Ilcweverj 
since tlie computa ciors are curried oijt perforce "’’ith rounded 
values, the relative ex’rors grow and overtake the i-ranted 
function when a straight forwcxrd use of the above recursion 
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iFj made in forward direction ( increasin,_, s) = i'hiB occurs 
vrien a lecurrence relrtion has t'‘o independent solu.tions 
and bhe solution de^.ix'ed is dij nish-iiig as the index 's’ 
increases, vliilo the compnaion sol'cbion is increasing. On 
reversinj' the dire c cion ':le ro] es of the buo solutions a-re 
interchuiigec and tcio contribution oi desired solution now 
increases while the un ’•anted so.liitjon Jim’ni.jhes ( Abramoifitz 
and otegun 1965). Compute. cion of the integrals Ip^(p) is 
th.'.^refore done applying the recur'sion in 'backward 

direction (decreasing s) . i’lie recursion is started from a 
sufficiently higher index ‘s’ using arbitre ,:-y starting 
values killer 1952) . The values obtained in this 

manner differ from the desired solution by a constant 
multiplier which is calculated from the vaJues of Iq(P) 
obtained from the tables for the error functions. 
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